Finite size emptiness formation probability of the XXZ 

spin chain at A = — ^ 



October 12, 2011 
Luigi Cantinill 



Abstract 

In this paper we compute the Emptiness Formation Probability of a (twisted-) 
periodic XXZ spin chain of finite length at A = — ^, thus proving the formulas 
conjectured by Razumov and Stroganov [3113 • The result is obtained by exploiting 
the fact that the ground state of the inhomogeneous XXZ spin chain at A = — ^ 
satisfies a set of qKZ equations associated to Uq{sl2)- 
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1 Introduction 



The computation of correlation functions in quantum integrable systems, is in general a 
quite hard task. One paradigmatic example is the spin-i XXZ spin chain pQ. Even in the 
study of the free fermion case quite interesting mathematical structures have appeared [2] . 
Starting from the seminal papers of Razumov and Stroganov [3],|1], we know that the anti- 
ferromagnetic ground state of the XXZ spin chain at the value A = — | of the anisotropy 
parameter (or equivalently q = e^'^*/'^) presents a remarkable combinatorial structure. The 
spin chain hamiltonians with an odd number of spins N = 2n+l and periodic boundary 
conditions or even number of spins N = 2n and twisted-periodic boundary conditions are 
related through a change of basis to the Markov matrix of a stochastic loop model [5]. In 
the loop basis, the "ground state" is actually the steady state probability of the stochastic 
loop model. The most astonishing discovery was made by Razumov and Stroganov [6], 
who observed that, once properly normalized, the components of the steady state are 
integer numbers enumerating Fully Packed Loop configurations on a square grid. This 
conjecture has been eventually proven in [7j. 

Despite being the XXZ spin chain at A = — | a fully interacting system, several of 
its correlation functions have simple exact formulate even at finite size. This is the case 
of the Emptiness Formation Probability (in short EFP), which is the probability that 
k consecutive spins are in the up direction in a chain of length A^. In their original 
papers [3lll], Razumov and Stroganov have conjectured exact factorized formulas for the 
EFP in terms of products of factorials. The aim of the present paper is to prove these 
conjectures. 

The XXZ spin chain with odd size has two ground states '^2n+i ^^'^ ^2n+i5 related by 
a spin flip on each site. Razumov and Stroganov have conjectured [3] that E2n+i{k), the 
EFP of a k-string of spins up in the state '^2n+iy satisfies 

- 1) (2A: - 2)\{2k - l)!(2n + k)\{n - k)\ 

E2n+i{k) {k - l)\{'ik - 2)\{2n - k + l)\{n + k - l)V ^' 

Strangely enough, Razumov and Stroganov didn't provide the analogous formula for the 
state '^2n+ii which reads 

EUijk-l) ^ (2fc-2)!(2fc-l)!(2n + fc)!(n-fc + l)! 
EL+i{k) {k-l)\{'ik-2)\{2n-k+l)\{n + k)\ ^' 

In particular, the probability of having a string of spins- up of length n (or n -|- 1) in a 
chain of length 2n -|- 1 is equal to the inverse of Anri'^n + 1), the number of of Half Turn 
Symmetric Alternating Sign Matrices of size 2?t, + 1, 

E2n^,{n) = El^,{n + 1) = A^^n + 1)"^ = f[ ^^^.".(g^^ (3) 

In the case of a spin chain with even length and twisted boundary conditions, the 
ground state is unique. Razumov and Stroganov have conjectured [1] that £'!„(/;;), the 
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Figure 1: Domain tiled by /c-Punctured Cyclically Symmetric Self-Complementary Plane 
Partitions, with an example of tiling for n = 3 and k = 2. The shadowed region indicates 
a fundamental domain. 

EFP of a k-string of spins up satisfies 

El^{k - 1) _ {2k - 2)\{2k - l)\{2n + k- l)\{n - k)\ 
El^{k) ~ {k - l)\{3k - 2)\{2n - k)\{n + 

The previous equation implies that in the case k = n 

El^{k) = AHT{2n)-\ 



k-l)\ 



(4) 



(5) 

Unlike the ground states of the odd size chains, whose components can be chosen to be 
real, the even size ground state has complex valued components, therefore we can consider 
also E2n{k), a sort of "pseudo" EFP obtained by sandwiching the ground state with itself 
(and not with its complex conjugate). The ratio of "pseudo" EFPs corresponding to the 
same size of the spin chain has a factorized form given by 



ELik) 



{2k - 2)\{2k - l)\{2n + k- l)!(n - k)\ 
{k - l)\{3k - 3)\{3k - l)(2n - k)\{n + k - l)V 



(6) 



It turns out that, apart for the factor — g, the ratio in eq.([n]) can be written as a ratio of 
enumerations of /c-Punctured Cyclically Symmetric Self-Complementary Plane Partitions 
(PCSSCPP) of size 2n, i.e. rhombus tilings of a regular hexagon of side length 2n, which 
are symmetric under a vr/S rotation around the center of the hexagon and with a star 
shaped frozen region of size fc, as exemplified in Figure [T] Calling these enumerations 
CSSCPP{2n, k) we have 



EL{k 



1) 



ELik) 



CSSCPP{2n, k-1] 
CSSCPP{2n,k) 



(7) 



For k = n one obtains 



ELin) 



-qYCSSCPP{2n) 
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where An is the number of Alternating Sign Matrices of size n 

and CSSCPP{2n) is the number of Cychcally Symmetric Self-Complementary Plane Par- 
titions in an hexagon of size 2n. The enumerations CSSCPP{2n, k) are easily computed 
by applying a result of Ciucu [10] concerning enumerations of dimer coverings of planar 
graphs with reflection symmetry. This will be explained briefly in Appendix [Bl 

Some partial results concerning the (pseudo)-norm, and the EFP have been obtained 
in the literature. In [TT], by cleverly exploiting the relation between a natural degenerate 
scalar product in the loop basis and the usual scalar product in the spin basis, Di Francesco 
and collaborators have proven eq.(l3]) and eq.®. 

In the large N limit — )■ oo, the EFP have been studied Maillet and collaborators [T2] . 
Of course the conjectures ([II [21 14]) must coincide and must give for the thermodynamic 
limit of the EFP 

„„g,(,)^f^V"nrO-i/3)rO + i/3) 

V 2 / r(j - i/2)r(j + 1/2) * ' 

This formula has been proven in [12] by specializing to A = — ^ a multiple integral repre- 
sentation for the correlation functions, which is valid for generic values of the anisotropy 
parameter A. 



The most effective technique which has allowed to compute (partial) sum of com- 
ponents in the loop or in the spin basis has been pushed forward by Di Francesco and 
Zinn- Justin [8] in the context of the periodic loop model with an even number of sites. 
They introduced spectral parameters in the model in such a way to preserve its integra- 
bility, the original model being recovered once the spectral parameters are set to 1. In 
this way the components of the ground state become homogeneous polynomials in the 
spectral parameters, which satisfy certain relation under exchange or specialization of 
the spectral parameters. As noticed first by Pasquier [13] and largely developed by Di 
Francesco and Zinn- Justin [T3] the exchange relations satisfied by the ground state are 
a special case very much studied quantum Knizhnik-Zamolodchikov 

equations (qKZ) [16]. In [T7] the authors have applied this idea to the XXZ spin chain 
with spectral parameters and have shown that the properly normalized ground state of the 
spin chain with periodic or twisted periodic boundary conditions satisfies a special case of 
the Uq{sl2) qKZ equations at level 1. Here we employ this property in order to compute 
the Emptiness Formation Probability. Our main idea is to consider a generalization of the 
EFP with spectral parameters EFP, which is constructed from the solution of the Uq{sl2) 
qKZ equations for generic q (see eqs. fl48ll49p . This "inhomogeneous" EFP has certain 
symmetry and recursion properties that completely fix it in the same spirit as the recur- 
sion relations of the 6-vertex model with Domain Wall Boundary Conditions completely 
determine its partition function. This will allow us to present an explicit determinantal 
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formula for the inhomogeneous EFP valid at g = e^'^^l'^ ^ and upon specialization of the 
spectral parameters will allow to obtain the formulas dH El HI [6]) . 

The idea to use the solution of the qKZ equation to compute the inhomogenous version 
of a correlation function can be in principle adapted to several other models like: XXZ 
spin chain with different boundary conditions, fused XXZ spin chain, Uq{sln) spin chain 
or even XYZ spin chain etc. [201423] . Indeed, in all these cases, by properly tuning the 
parameters (generalizing the relation q = e^'^*/^), one has a so called "combinatorial point" 
at which the ground state energy per site doesn't get finite size corrections. By reasonings 
similar to the one in [T7], one can argue that the ground state with spectral parameters 
satisfies a qKZ equation (or, in the case of the XYZ spin chain, an elliptic version of it). 

Whether this idea could lead to other exact finite size formulae for some correlation 
function is an open question that in our opinion deserves further investigation. 

The paper is organized as follows. In Section [21 after having recalled some basic facts 
about the XXZ spin chain, following [T7] we present the exchange equations satisfied by 
the ground state at A = — |, then in Section we derive the recursion relations satisfied 
by the solutions of the Uq{sl2) qKZ equations at level 1. In Section [3] we define the 
inhomogeneous version of the (pseudo) EFP, constructed using the solutions of the qKZ 
equations. We derive first its symmetries and then in Section 13.11 we derive the recursion 
relations which completely determine it. In section S] we will restrict to g = e^'^*/^ and, 
by showing that certain determinantal expressions satisfy the same recursion relations 
as the inhomogeneous EFP, we produce a representation of this inhomogeneous EFP 
whose homogeneous specialization is considered in Section 14. ![ where we prove the main 
conjectures. In Appendix [A] we compute the determinants of a family of matrices which 
are relevant for the computation of the homogeneous specialization considered in Section 
14.11 In Appendix [B] we compute the lozenge tilings enumerations CSSCPP{2n, k). 



2 XXZ spin chain at A = — | and qKZ equations 

The hamiltonian of the XXZ spin chain acts on a vector space Hn = (C^)*^^ that consists 
of copies of each one labeled by an index i. The hamiltonian is written in terms of 
operators af which are Pauli matrices acting locally on the i-th component Cf 

hn{^) = "2 E + ^r<i + (11) 

i=l 

It is convenient to parametrize the anisotropy parameter as A = — . The model is 
completely specified once the boundary conditions are provided. Here we will consider: 

• periodic boundary conditions for odd values of the length of the spin chain, A^ = 
2n + 1, i.e. a%^^ = tr" 

• twisted periodic boundary conditions for even values of the length of the spin chain, 
A^ = 2n, i.e. cr^_^_i = cr^, while (t^_^_i = e^^^a'^j^^, where a"^ = + ±.ia'^ . 
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It is well know pLj that the Hamiltonian ( TTT]) . for generic values of the parameter A and 
of the twisting, is the logarithmic derivative of an integrable transfer matrix. In order 
to define the transfer matrix we need the i?-matrix and the twist matrix. In the present 
context the i?-matrix is an operator depending on a spectral parameter which acts on 
a tensor product Cf (g) Introducing the basis of 



we write Rij{z) in the basis {ej ® ej, e\ 









with 



a{z) 



qz — q 
q — q~^z 



b{x) 



z-1 

— q '^z 





biz) 

C2{z) 







1 

) e}, e| 


ci{z) 
b{z) 




8 ej} of Cf 

\ 




a{z) j 



C] as 



(12) 



ci{z) 



q — q~^z 



C2(Z 



{q-q 1) 



q-q- 



The twist matrix Q{4>) acts on a single and in the basis (ej, e^) it reads 



^(0) 



(13) 



(14) 



Using both the twist and the i?-matrix we construct the family of transfer matrices 

7V(l/|z{i,...,Ar},0) = tro [Ro,iiy/zi)Ro,2{y/z2) ■ . ./2o,iv(y/^jv)^o(0)] (15) 

depending on "vertical" spectral parameters Z{i 7v}0 and a single "horizontal" spectral 
parameter y. Thanks to the commutation relation 

[/2,,,(x),a(0)®l^,(0)] = 

and the Yang-Baxter equation 

Ri,j{zi/ Zj)Ri^k{zi/ ZkjRj^Zj/ Zk) = Rj^Zj / Zk)Ri^j^{zi/ Zk)Rij{zi/ zj) 
the transfer matrices for different values of y commute 



[TN{y'\ 



'{1,...,N}, ' 



^Jv(y'1z{l,...,Ar}, 



0. 



The hamiltonian of the XXZ spin chain is given by 



rfTjv(2/|l, 



T^(l|l,' 



dy 



y=l 



q-q- 



3N 

HNiA) - —A 



(16) 
(17) 
(18) 

(19) 



At A = — I and for generic values of the vertical spectral parameters, both in the odd size 
case with periodic boundary conditions and in the even size case with twisted boundary 
conditions, the transfer matrix has an eigenvalue equal to 1. 



"'^Our convention for a ordered string of variables labeled by an index is to use a bold character and a 
label for the ordered set of indices of the variables: ^{ai,...,aN} — {^ai j ■ • ■ j ^a^}- Often, when clear from 
the context, we will omit the label {ai, . . . , a^} and write x for X{qj^_ (jj^j. 
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• When = 2n+l, the eigenspace with eigenvalue 1 is two-fold degenerate, \l/^„_,_i(z), 
corresponding to two the values ±| of the total spin = ^ J2iLi '^i 

5^vl/±^i(z) = ±vl/±^i(z) (20) 

The two eigenstates are related by a flipping of all the spins 

N 

*^nH-l(z)=n<-*2~n+l(z). (21) 

i=l 

• When N = 2n there is a single vector \E'2n.(^) with eigenvalue 1. It is in the zero 
sector of the total spin, 

S^n^iz) = 0, (22) 

These eigenstates reduce to the anti-ferromagnetic ground state (s) of the XXZ spin 
chain when the spectral parameters are specialized at 2;^ = 1. 

2.1 Exchange relations at A = — ^ 

A crucial observation made in [T7| was that, for an appropriate choice of the normalization 
of ^^(z), the eigenvector equation 

T^(y|l,0)vl>^(z) = vl>^(z) (23) 

is equivalent to a set of exchange relations. Define the exchange operator -Pij (ef ® ej) = 
<S)ej, the left rotation operator (t{vi ® f 2 ® ■ ■ ■ ® v^-i ® f at) = V2^ - ■ ■ ^v^^i ® f at ® f 1 
and let Ri^i^i{z) = Pi^i^iRi^i^i{z) , then \E'^(z), as a polynomial of minimal degree in 
the spectral parameters Zi, is determined up to a constant factor, by the following set of 
equations [17] 

Ri,i+i{zi+i/ Zi)'^'^{zi, . . . , Zi, . . . , 2; at) = \E'^(zi, . . . , Zi, . . . , zn) (24) 

(T^^(Z1, Z2,..., Zn-1, Zn) = D'^%{Z2, Zn-1, Zn, SZi). (25) 

with D = s = 1. These equations can be seen as the special case q = e^'^'^l'^ of the level 1 
qKZ equations [16], which corresponds to generic q, s = (f and D = g3A^g3(s^+i)/2_ rpj^^ 
solution of the level 1 qKZ equations can be normalized in such a way that they become 
polynomials in the variables Zi [T7| of degree n — 1 in the case of even size = 2n, and 
degree n in the odd case A^ = 2?7, + 1. 

Using the projectors pf = of the i-th spin in the up/down direction, let us write 
the exchange eqs. (124|) in components. If we have a spin up at site i and a spin down at 
site i + 1 or viceversa, we have 

+ j f,, . _ + {qzj - q-hi+i)'^''{zi+i, Zi) - (g - q-^)zi<il''{zi, Zi+i) 
Pi Pi+i^ [Zi,Zi+i) — cr. (J-^^ — 

(26) 

+ - N + - {qZi-q-hi+l)'^^'{zi+l,Zi) -{q-q^^)zi+l^^'{zi,Zi+l) 

ptp,+i'^^{z^,Zi+i) = (T+(T,.+i 

Zi+l — Zi 



These equations form a triangular system. Starting from a given component we can 
reconstruct all the others repeatedly using eqs. fl26l) . Therefore if we want to show that 
two a priori distinct solutions of the qKZ equations actually coincide, it is enough to check 
the equality of one of their components. 

When there are two consecutive spins pointing in the same direction at positions i and 
i + 1, then the first of the qKZ equations reads 

pfpt,<i^>'i. . . , ^.+1, . . . ) = '^''^''\''' pfpti'^^{. . . , ^.+1, . . . ) (27) 

qzi- q ^Zi+i 

which means that the components having two consecutive spins up or down at positions 
i and z + 1, pfpfj^^'^^{. . . , Zj+i, Zi, . . .) have a factor qZi — q~^Zi+i 

pfpf+i'^^ ■■,z,, . . . ) = {qz, - ..,Zi, z^+i, . . . ) (28) 

and the vectors ■ ■ ,Zi, Zi+i, . . . ) are symmetric under exchange Zi -H- Zi^i. Another 

useful relation is obtained by considering the matrix Cj, which is proportional to a projector 
and is a generator of the Temperley-Lieb algebra 



/ \ e? = re„ t = -q - q'^ 

-g 1 

1 -g-^ ' 

yo / eiej=ejei for |z-j|>l 



eiCiiiCi = Ci (29) 



The matrix Cj is preserved under multiplication by a i?-matrix for any value of the spectral 
parameter 

eiRi^i+i{z) = Ri^i+i{z)ei = e^. (30) 
By applying ej to the left of the first of the qKZ eqs. fl24p we find 

ei^f^(. . . , 2;,, . . . ) = ei^"^(. . . , z,+i, Zi,...). (31) 

The components with most consecutive aligned spins have a completely factorized form 

j e / N _ TT g^i - q'-^Zj -p-r qzj - q~^Zj 

^t,...A,l,...,i^^^- 11 g-g-i 11 g-g-i 

V V ' l<i<j<n ^ ^ n+l<i<j<2n ^ ^ 

n n 

_1 _1 2n+l 

,Tr+ ^ N TT ~ 1 TT ~ 1 TT 

*t,...,u,...,i^^)= n -j^^ n ^3^n-^ (32) 



n + l 



l<«<i<»i+l n+2<i<i<2n+l ^ i=n+2 

qZi - q ^Zj l^i ~ 



l<i<j<n ^ ^ n+l<i<j<2n+l ^ ^ 



n+l 



where the residual normalization ambiguity has been fixed by requiring these components 
to be equal to 1 for Zi = 1. From eqs. (l32|) we see that the maximally factorized components 
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satisfy (among others) the following relations 

J ; • • • ; h^^; ' ' ' ^ J ^ ' ^ ' ^' ^ 

(33) 

n |(z) = (l-g')" lim ^1 |(z) 

n n „ n+1 

Using the triangularity of the eqs. we can conclude that the eqs. ([32]) induce equalities 
between components of '^2n+ii^) ^2n(z) and components of '^2n+2i'^) ^2n+i(z) with 
the last spin down, i.e. 

^2'n+l(z) ® 4„+2 = (1 - '?~')>2~n+2^2n+2(z)|.2„+2=0 

*2n(z) ® 4n+l = (1 - 22"n+lP2"n+l*2"n+l(z) ^^^^ 

We will use these equations in order to find a relation among the inhomogeneous versions 
of the EFP that we shall introduce in Section El 



2.2 Recursion relation 

We claim that, upon specialization Zj+i = g^Zj the solution of the qKZ equation for 
spins reduces to the solution of the same system of equations for — 2 spins. In order to 
make the previous statement more precise we need to introduce some notation. Let Vi be 
the vectors which are in the image of the projectors proportional to the generator of the 
Temperley-Lieb algebra Cj 

= ej®e|+i -g"^e|®ej+i, = -{q + q~^)vi (35) 

Introduce the injective map 

$» ; (c2)«^ ^ (C2)®^+^ (36) 

which inserts the vector Vi at position + 1) and shift by two steps the indices of the 
sites with j > i, i.e. on a basis 

^^SieT ® ■ ■ ■ ® e^ri' ® ® ■ ■ ■ ® e^^) = 6"^^ ® ■ ■ ■ ® e^'i ® Vi ® e'^^^ ® ■ ■ ■ e^^^- (37) 
Then we claim that 

Proposition 1. The solutions of the exchange equations [24\ ), with the "boundary con- 
ditions" given by eqs. satisfy the following recursion relations 

i-l _i 3 _i 

j=l ^ ^ j=i+l ^ ^ 

_ (38) 

with 5(e) =6{-) = 0, 6{+) = 1, f^ii) = ^ - Lf J = i^i' • • • > ^i+iy • • • > ^n}, 

i.e. the ordered set z from which the variables Zi and Zi^i are removed. 
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Proof. If Zi+i = q^Zi then Ri^i^i{zi/ ZiJ^i) becomes a projector proportional to a generator 
of the Temperley-Lieb algebra 

RiM^l''') = ^"'e, (39) 
Therefore, by specializing the qKZ equations to Zj+i = g^Zj, we deduce 

-^^j^{zi, Zi+i = q^Zi) = Ri^i+i{q~'^)'^%{q'^Zi, Zi) = r'^d'^'^iq'^Zi, Zi). 

In particular '^'^{zi, Zi+i = q^Zi) lies in the image of $j^''_2 and, by the injectivity of this 
map, there is a unique \&^(2;j, Z||--pjj) such that 

v|/i^(^„ ^,+1 = q^z^) = (40) 

In order to determine the equations satisfied by '^^^^{zi^Ti^^rjj-^^) we make use of the fol- 
lowing relations among i?-matrices 

ei^i_i(zi+2/^i)A(^i+2/(g^^i))^i+i(^i+2/^i-i)^i(g^^i/2;i-i)Aj„i(zi/2;i)ei = 

- q~^Zi){qZi - q'^Zj^i) ^ (4^) 
[qZi^i - q~^Zi)[qZi - q-^Zi+2) 

Applying both sides to '^'^{zi, ^,+1 = q'^Zi) and using that 

Ri-l^i+2{Zi+2, Zi)^i = ^iRi^i^i{Zi+2, Zi) (42) 

we get that ^^^(2;,, z^r— pj^) satisfies 

{qZi+2-q~^z,){qZi-q-^Zi_i) g, ^ ^ ^ ^^^^ 

[qzi_i-q ^Zi){qZi-q ^Zi+2) ^'^^ 

and the vector 

(-g)-/-(-)(g^.0^^^^n '"'1 n 3'"' 1 (44) 
^^qzi — q '-Zi q-^Zi — q ^Zj 
j=i j=i+i 

satisfies all the qKZ equations at size N — 2. In order to check that it coincides with 
\l/^_2(Z|j-r3pjj^) it is enough to check that the components with most ahgned consecutive 
spins starting from position z + 1 coincide, which is indeed the case. □ 



3 (Pseudo)-EFP with spectral parameters 

The formal definition of the Emptiness Formation Probability makes use of the natural 
scalar product (■, ■)n on Hat, induced by the scalar product on where {ej, e|} form an 
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orthonormal basi^. The (pseudo)-EFPs read 

fe+i(i),nliR^-^rn+i(i)) 

ELik) 



(^2n(i),nliR^-^ 



(45) 



(vl/|„(l),vl/^„(l)) 

Our strategy to compute the EFPs is to consider an inhomogeneous version of these 
quantities which is obtained, roughly speaking, by substituting in eqs.( HH]) with 
\E'^(z), solution of the qKZ equations. We shall see that if the substitution is done in the 
proper way, then the inhomogeneous EFPs turn out to be symmetric polynomials in the 
spectral parameters and satisfy certain recursion relations which completely characterize 
these functions among the polynomials of the same degree in the variable Zi. 

When defining the inhomogeneous EFP for k aligned spins up, it is convenient to 
extract the factor Ui^i^j^kilVi - Q'^Vj) from IlLi Pi^^A^(y{i,-,fc}; '^{i,...,N-k}), and to in- 
troduce the vectors \E'Ar(/i;; y{i,...,fc}; Z{i,...,Ar_fc}) e HN-k 

f, X nti Pt'^Njyi, ...,yk,zi,..., ZN-k) 

^" "-"^ = — n.,.<,,.ta-rt) — 

Let us moreover introduce the operator 

N 

VN{z) = l[{z,pt+p;). (47) 

1=1 

that multiplies each component of the vector by Zi for a spin-up at position i. The last 
ingredient we need is the * operation which consists in substituting q with q~^. Our 
definition of the inhomogeneous (and unnormalized) EFP is the following 




^N-k ,-^(A^)(p^_^,(,)(vi;^(fc; g"Vu.+i,...,2.}; z))*, Kik; y{i,...,.}; z)) 



^Nik',y{i,...,2ky,'Z{i,...,N~k}) = 

(4J 



N-k 



where S{e) = 6{—) = 0, while 5(+) = 1 and has the parity corresponding to fi. For 
the moment it is evident that for /i = e or /i = — , S'^{k;y;z) are polynomials in their 
variables. Actually the polynomiality is true also for the case // = + as will be shown 
below. We define also the inhomogeneous version of the pseudo-EFP 

£Nik;y;z) = 

n-=.+i 1/!"^^"' nf=f zl'^^-'iKik; g- V+i,...,2.}; Kik; y{i,...,.}; z))^.,. 

(49) 



^In the following most of the time it will be clear from the context which Hilbert space we are 
considering and therefore we will omit the label N in the scalar product (•, •)jv- 
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The choice to multiply the variables i/k+i, ■ ■ ■ ,y2k by is motivated by the fact 
that in this way £'^{k;-y;z) turns out to be symmetric under exchange i/i -H- i/j for all 
1 ^ "^j J < '2k, as will be shown at the end of Section [STTl The polynomials (148| |49|) have 
other remarkable properties, but for the moment let us simply notice that for q = e^'^*/^ 
and Zi = yi = 1 these functions reduce to the unnormalized version of the EFPs as defined 
in eqs. (H5|) . 



Apparently, looking at eqs. fl48f49p . it seems that we have six different families of polyno- 
mials under consideration. Actually in the odd size case we have 

^2n+iik; y; z) = S^^^^{k; y; z). (50) 

This follows from the fact that 

2n+l 

n =^2n-M(z)(^2V(^))* 

(51) 

2n+l ^ ^ 

n ^r*2n+l(z-') =7^2n+l(z)(vi/2-+i(z))*. 

i=l 

To prove eqs. flSTj) we observe that the vector \E'^(z~^) satisfies the exchange equation 

z~') = i?,,H_i(z,/z,+i)vl/^(z-\ z-\). (52) 
The same equation equation holds also for the vector 'P7v(z)(\E'^(z))*, i.e. 

VN{zi+i, Zi){'if'^{zi+i, Zi))* = Ri^i+i{zi/zi+i)VN{zi, Zi+i){'^'^{zi, Zi+i))*. (53) 

This is a consequence of the following commutation relation among the i?-matrix and the 
operator {pf + Zip-){pt+i + ^i+ip'+i) 

Ri,i+i{zi/zi+^){zipt +p~){zi+ipt+^+p^^^) = {zi+ipf + p-){zipt^^ + p-^^)R*i^^{zi+i/ Zi). 
which implies 

Ri,i+i{zi/zi+i)V{zi, Zi+i) = V{zi+i, Zi)Rli^^{zi+i/ Zi) (54) 

and eq. (l53|) . Therefore to conclude eqs. llSTl) it is sufficient to check that they hold for the 
components with most aligned spins. 

We are left with only four different inhomogeneous EFP and as a bonus we have also 
shown that i^2ri+i(^j ^) ^ polynomial of its variables. 

Symmetry under Zi O Zj 

The inhomogeneous (pseudo)-EFP £'^{k;y;z) is obviously symmetric under exchange 
yi O yj for 1 < i, j < k and k + 1 < i, j < 2k. Using eqs. fl52|53l) it is easy to show that it 
is symmetric also under exchange Zi O Zj. Indeed 

^Tvl^j y; • • • ) Zi+i, ■ ■ ■) = 

{(V{zi, Zi+i)'^%{k; Zi, Zi+i))*, Ri^i+i{zi/ Zi+i)^^(/c; Zi+i, Zi)^ = 
{Ri^i+i{zi/zi+i){V{zi, Zj+i)^^(/c; Zi, Zi+i))*, ^'^^(/c; Zi+i, Zi)^ = (55) 
V{zi+i, Zi)<i''^{k; Zi+i, Zi))*, ^'^(fc; Zi+i, Zi)^ = 

^nQ^'i y; • • • 5 ZiJ^i, Zi, . . .) 
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where in the third equahty we have used the fact that the ^-matrix is symmetric while 
the fourth equahty follows from eq. (l53|) . The proof of the symmetry of the pseudo-EFP 
under Zi O Zj is completely analogous. 

Factorized cases 

Using eqs. fl5^ we can provide the value of S'^^^{k;y;z) corresponding to the maximal 
number of consecutive aligned spins. They coincide for the true and for the pseudo EFP 
and read 



^2k['^^y^^)- 11 (a-a-^y^ 

l<i<j<k ^ ' 



l<i<j<fc ^ ' 1=1 

c~/-(l,.^.^^- rr {qzi-q-^z,){qzj-q-^z{) 
^2k+i\'^^y^^) - 11 (a_a~i\2 
i<i<j<fe+i ^ ' 

We will see in the following that the first and the third of these equations will provide 
the starting point of a recursion which will be worked out in the next section and which 
completely characterize the inhomogeneous EFP. 

3.1 Recursion relation for the inhomogeneous EFP 

We begin this section by presenting some relations among the EFP at different parities 
which are obtained by setting one of the spectral parameters to zero or sending it to 
infinity 

(2n+l-fe \ 
n ^rM^2V2(^;y;z)l..„+.=o (57) 

£',Sk- y; z \ = - q-'f" hm z^^l,E^^,^,{k- y; z). (58) 

The first of these equation follows from eqs.f l34p and by noticing that 

^(z) l^2n+2=0 = ^(Z \ ^2n+2)P2n+2- (59) 

For the second one we notice that, writing 

^2n+i(^;y;z) = 

((^2"n+l(^; Z))*, P2„,(Z \ Z2n+l)z2n+lVtn+l-^2n+l{k-, ^)) 2n+l + (60) 
((* 

2n+l(^jz)) ,^271(2 \ 22n+l)P2n.+1^2n.+l(^j z))^^^^ 

the first term in the r.h.s. is a polynomial in Z2n+i of degree 2n — 1 while the second is 
of degree 2n, therefore in the limit only the second one survives and we can use again 
eqs.( 



12 



Specialization Zi = q^'^Zj 

The inhomogeneous EFP satisfies a recursion relation inherited from the the recursion 
relation among solutions of the qKZ equations, eq. (l38|) 

^Af(^j yj z;)|z,_^i=g2^. = 
T~^{{^!^^{k]Zi,Zi+i = q~^Zi)y,V{zi+i = q^Zi)ei^!^^{k]Zi,Zi+i = q^Zi) = (61) 
T~^{eiV{zi+i = q^Zi){^!%{k] Zi, Zi+i = q~'^Zi)y, '^'^{k; Zi, Zi+i = q^Zi). 

A simple computation shows that ei{pf + Zip^) {pf +q'^ Zip^) = ^ei{pf + ZiP^) {pf +q'^ ZiP~)e* 
which means ^ 

eiV{zi+i = q^Zi) = -eiV{zi+i = q^Zi)e* (62) 
r 

and we can substitute it into the last line of eq. fl6ip obtaining 

T~^{eiV{zi+i = q^Zi){ei^!^j^{k] Zi, Zi+x = q'^Zi))*, ^'^(/c; Zi, Zi+i = q^Zi)N (63) 

Now use eq. fl3ip in order to exchange the variables Zi and ^j+i in the l.h.s. of the scalar 
product 

T~'^{eiV{zi+i = q^ Zi){ei^%{k; Zi+i = q~'^Zi, Zi))*, ^^(/c; Zi, Zi+i = q^Zi)^. (64) 



Therefore we can apply to both sides of the scalar product the recursion relation fl38l) and 
find ^ 



E'^_2{k;y;z\ {2:^,2:^+1}) 



2k 



,2^ T-fm^ZJL^ TT (g^j - g ^z,){q^Zi-q ^z^) 

j=l ^ ^ i<j<N-k ^ ' 

The case of the pseudo EFP at even size is analogous 



(65) 



^2n(^;y;z)|^,+i= 

(*^„(/c; z-^^ = q-^z-^), ei%^{k; z,, Zi+i = q^z{)2n = (66) 

and again we can apply the recursion at the level of vectors to both sides of the scalar 
product finding 

^In(^;y;z)l..+i=g2^, 



^2n-2i.k]y]2.\{Zi,Zi+i}) 



(-1)^(1+^^)17^^^^ TT ^^^^ 

Let us look at £^2n^(A;; y; z) as polynomials in Zi. Their degrees are in both cases less 

than 2n — 1. The recursion relations eqs. fl65|67l) provide the value of £^2n^(/c; y; z) for 
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2(2n — k — 1) distinct values of zi (i.e. for zi = q^Zi). Therefore, for n > /c, by Lagrange 
interpolation these specializations determine uniquely £^2n^(A;; y; z) once £^2n-2(^! ^) 
known. 

As a first consequence we can argue that £2n{k; y ; z) is symmetric under exchange yi ^ 
yj for all 1 < < 2k. Indeed for the case when n = k we have explicit expressions for 

^2n=2fc(^j yj given by eqs. (l56i) . from which we can read that they are even independent 
from y. The recursion relations fl65f67p are symmetric under exchange yi -H- yj and 
therefore by induction, if S2n_2{k; y; z) is symmetric, then also S2n{k;y;z) is symmetric. 

A second important consequence is that any family of polynomials labeled by n and 
k, which satisfy the following conditions: 

- they are symmetric in the spectral parameters, 

- the degree in each spectral parameter is less than 2n — 1, 

- they coincide with 8l^{k\ y; z) for n = k, 

- they satisfy the recursion relations eqs. (16511671) . 

must coincide with S2n{k; y; z). This line of reasoning will be adopted in Section H] where 
we will provide a determinantal representation of £2n{k\ y; z) at g = e^'^*/^. 

The same arguments holds also for £^2n+i(^j '^)^ because the degree is less than 2n + l 
and we have always enough specialization in order to apply the Lagrange interpolation 
and reconstruct all the £^2^^_^(A;; y; z) starting from the initial conditions £^2fc+i(^i yi 
The case of S^n+A^'i 2) is slightly different. Again the degree is bounded by 2n + 1 and 
this allows to fix £^2^1+1 (^i Yi z) starting from £^2^+1 (^j Yi z), but the problem is that we 
do not have an explicit formula £^2A:+i(^) Yi 2:), being available only ioi n = k — 1. This 
apparent problem is bypassed using relations (IFHl . 



4 Inhomogeneous EFP at q-^ttz/s 

In order to introduce the expression of S'^{k; y; z) and of £'{^{k\ y; z) which is best suited 
for taking the specialization Zi = y^ = ^ we analyze first the case /c = 0, in which there 
are no variables y. Let us introduce the Young tableaux 

r,r , ,r + 1 , ,r + z — 1, ,r + m— 1,^ , , 

A(m, r) = { L-J , L^J , . . . , L ^ J , • • • , L ^ J } (68) 

Then we find that the inhomogeneous version of the squared norm or of the sum of the 
square of the components is given in terms of the product of two Schur polynomials 

£j^(0; z) = 3-Lf J(rf . . . , ZN)SxiN,i){z,, . . . , ^^^^ 

^IJO; Z) = 3-"("-l)5A(2n,0)(^l, • • • , Z2n? 

The proof of eqs. fl69p is quite simple and follows the pattern discussed at the end of 
the previous section. Eqs. (l69|) are trivially true for = 1,2 (or n = 1), moreover 
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their r.h.s. are polynomials in zi of degree at most 2\y~\ — 1- The Schur polynomials 
Sx{m,r){zi, . . . , Zm) Satisfy a recursion relation when one specializes Zi = q^Zj (see for 
example Appendix B of [19]) 



X{m,r) {'^)\zi=q±Zj 



-(f^^iY n ~ ^"^^j) 'S'A{m-2,r)) (Z \ {Zi, Zj}). 



(70) 



1=1 



This implies that the r.h.s. of eqs. (l69|) satisfy the recursion relations eqs.( l65|67|) and 
therefore eqs. (l69|) hold. 

Generic value of k 

The recursion relations eq.f lTD]) for the Schur functions Sx(m,r) suggests a possible repre- 
sentation also in the case k ^ 0. For the sake of clarity let us focus for a moment on 
£|„(/c; y;z). It is easy to see that any product of the kind S'A(2n,o)(y/5 z)S'A(2n,i)(y/'=, z), 
with I C {1, . . . , 2k} and I'^ = {1, . . . , 2k} \ I, satisfies the recursion relations (165|) . but 
with a "wrong" initial condition. It is reasonable to hope that an appropriate linear com- 
bination of terms with different choices of / could provide the right initial condition and 
hence £|„(/c;y;z). 

In order to present how this idea actually works it is convenient to introduce a bit of 
notation. Let p, a be strictly increasing infinite sequences of non negative integers, then 
consider the following family of matrices 

/ 



4' 




■ l 





. 







zl^ . 


• ^2 





. 














. 








. 





4' 


4' ■ 







. 





4' 










. 



















y? ■ 


■ y'r- 


yl' 


yl' ■ 


. yl'^' 


y'2' 


y? ■ 


. y^' 


\ yl\ 


yt ■ 


■ vt' 


yl\ 


y?s ■ 


• y2s 



(71) 



and let us define the following polynomials 

det A^('^''^)(r,s;y;z) 



5(P,<^)| 



r,s;y;zj 



ni<j<i<rV 



Hi 



<i<3<2s 



iVi - Vj) n i<i<r (Zt - Vj)' 



(72) 



The divisibility of det M'^P'^\r, s; y; z) by Hi 



<i<j<r \ 



and by Hi 



<i<j<2s 



iVi - Vj) is 



immediate. If we set Zi = yj for some i, j then we subtract from the row corresponding to yj 
the two rows corresponding to Zi getting a null row. This means that det A^*-^''^^(r, s; y; z) 
is also divisible by {zi — yj). 



15 



Using the Laplace expansion along the first r + s columns we can write S^^'^'^r, s; y; z) 
as a bilinear in Schur polynomials 

n n (yi~yj) 

S'-P^ '{r,s;y;z) = ^ (-1) ^ ' Sp(^+,)(z, y7)5^(^+,)(z, y/c 

/c{i,...2s} 11 KV^-Vo) 

l^l=s l<i<j<2s 

(73) 

where pirn) and (T(m) are Young tableaux of length m, whose entries are p{m)i = pi — i + 1, 
a{m)i = cTj — i + 1. In particular notice that when s = then iS*-^'°'^(r, 0; z) factorizes as 
product of two Schur polynomials. 

Now let us introduce the following family of integer sequences 

A(0) = {0,1,3,4,6,7,...} 

~X(r)-\^l^l±l\ A(l) = {0,2,3,5,6,8,...} 

2 J' A(2) = {1,2,4,5,7,8,...} ^'"^ 

Then we claim that 

S^n^k- y; z) = 3-"^+'=(^-i)/25(^W'^(^))(2n + 1 - fc, A:; y; z) (75) 

(2n-fe+l \ 
n ^7^j'5(^^^^'^^^))(2r. + l-A;,A:;y;z) (76) 

S^^ik; y ; z) = 3-n(n-i)+fc(fe-i)/2^{A(o),A{i)) - k, k; y; z) (77) 

(2n-fc \ 
JJ ^-1 5(A{0),A(2)) - A;, A;; y ; z) (78) 

These formulas reduce to eqs. fl69l) for /c = 0. Using the relations among inhomogeneous 
EFP with different parities eqs. fl57f58p and the explicit form of the matrices we 
see easily that eq. (!76|) follows from eq. (!77|) , which in turn follows from eq. (175|) . Therefore 
it remains to prove only eq. (17^ and eq. (l75il . 

The r.h.s. of eq. (|75]) and eq. (175i) are polynomials in Zi respectively of degree 2n + 1 
and 2n — 2. Moreover using eq. (1701) and the form of ^'•'^'•^^•'''^'•^^^^(m, k; y; z) expressed by 
eq. (l73p we can easily obtain the recursion relation 

{-q^z,r+^^ n - ^^^^■)' n(^- - g^^.)5(^('-^)'^(^^»("^ - 2, A;; y; z \ {z,, z,}). ^^^^ 

i=l a=l 

In order to conclude, as explained at the end of Section 13.11 it remains to show that 
eqs. (l75ll78p hold for n = k, i.e. we need to prove that 

5(A(i).A(2))(^ + 1, A;; y; z) = ni<.<,<.-,i(^.^ + + ^]) ...^ 
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We proceed by factor exhaustion. A preliminary remark is that both S^^^^^'^^'^'^\2n — k + 
l,A;;y;z) and S^^^^^^^^^'^\2n - k,k;y;z), as polynomials in yi are of degree n — k and in 
particular they vanish as soon as /c > n. Therefore using the recursion relation fl79|) we 
conclude that 

5(A{l),A(2))(^^l^^.y.^)|^^^^^^^^5(W^^ (81) 

Since their degree as polynomials in Zi is respectively k and /c — 1, this means that we 
have proven eqs. llHT]) up to a numerical constant. Such a constant will be fixed to be equal 
to 1 in the following section, where we shall compute explicitly the specialization of the 
inhomogeneous EFP for Zi = and yj = t^-^-'+i-i. 



4.1 Homogeneous limit 



In this section we arrive at last to the computation of the homogeneous (pseudo) EFP 
using eqs. fl75ti75]l . We need only to consider a last intermediate step by setting z = z(t) 
and y = t^~'^y(t) with 



and y(t)j = 



The matrices 7\/l(A(r),A(s))(^^ ^. 2,{ty^ t'"y(t)) with r, s and m as in eqs 
able structure as columns matrices. Let us look at a concrete example 



have notice- 



_y^(A(o),A(i))(3^2;z(t);tVW) 



/ ^0-0 


^1.0 


^3-0 


^4.0 


^6-0 

















\ 


^0.1 


ill 




^4.1 


^6.1 



















^0.2 




^3-2 


^4.2 


^6-2 


































^0-0 


^2.0 


^3-0 


^5-0 


^6-0 



















^0.1 








^6-1 



















^0.2 


^2.2 


^3-2 


^5-2 


^6-2 




^0.3 




^3-3 


^4.3 


^6-3 


^0-3 


^2.3 


^3-3 


^5-3 


^6-3 




^0-4 


^1.4 


^3.4 


^4.4 


^6-4 


^0-4 


^2.4 


^3-4 


^5-4 


^6-4 




^0-5 




^3-5 


^4-5 


^6-5 


^0-5 


^2-5 


^3-5 


^5-5 


^6-5 




1^ 


^1.6 


^3-6 


^4-6 


^6-6 


^0-6 


^2-6 


^3-6 


^5-6 


^6-6 


/ 



^2) 

The entries of the j'-th column (apart for the zeros) are consecutive powers of some Vj, 
where Vj is itself a power of t which depends on the column index j. In the precedent 
example v = {1, t, t^, 1, t*'}. Moreover some Vj appear twice, once in the 

first half of the columns and once in the second half (in the example l,t^,t^), while the 
remaining Vj are of the form aiA* in the first half of the columns and 02 A* in the second 
half (in the example A = t^, ai = t, 02 = t^). 

By these considerations we are led to introduce the following families of 2(£ + r + s) x 
2(£ + r + s) matrices ^(^'''''^^(v; A, 01,02), that are made of 6 blocks of rectangular matrices 
as follows 



'(£,r,s) 



v; A, 01,02) 







D 



£+r;r+s('^l^) 





D 



[aiX] D 





(0) / 
2s-l \ 



D 



(0) 



\ 

02A) 

(02A) J 
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where v = {vi, . . . , vg}, ai\ = {oj, OjA, . . . , a^A''"'"'^ ^} and the blocks consist of the follow- 
ing rectangular matrices 



( v{ ... "^^ \ 

j + 1 7 + 1 7 + 1 

V\ V2 . . . 



V ^1 ... vi ) 



^4) 



^5) 



Apart for a trivial reordering of the columns we have 

_;^(A(0),A(l))^2n - k, k; Z(t); y(t)) = gin,n-k,k)^^^3i~3y^ ^3^ ^2) 
_y^(A(0),A(2))(2^ - fc, /t; Z(t); y(t)) = ^(->n-fe,fc)({t3.-2|. ^3^ ^2) 

;U(A(o),A(i))^2n + 1 - A;,A:;z(t);y(t)) = ^("+i'"-'='^)({t3i-3}. ^3^ ^2^ ' 

_;^(A(1),A(2))(2^ + 1 - fc, A;; z(t); y(t)) = gin,n-k+l,k)^^^3i~2y^ ^3^ ^) 

Therefore, by calling 

SUk;t):=£^{k;y{t);zit)) (86) 

we have 



Ce .^..X _ ^-n(n-l)+k(k-l)/2 det(g("--'-'')({i3»-3};t3,M^)) 
Ce _ o-n(n-l)+fc(fc-l)/2 det (g(->"-'^.'=) ({t3--^ };i3,l,t2)) 

<^2nl'^i ''J — -J n , to-i-t'-im .,(P 



(87) 



ri-- - 0-n^+k(k-l)/2 det(g("+i>"-''-'')({t3^-3};t3,t,f2)) 

<^2n+ii'-,^j - ^ ni<,<,<2„-fe+i(*^-^-«'-^)ni<.<,<2n+fe+i(t^-^-t'-^) 

The remarkable fact about the matrices ^'-^'^'^■'(v; A, ai, 02) is that their determinants 
factorize nicely. For r > we have 

det6;(^'^'^)(v;A,ai,a2)= J] " ^j)' H H (v^ - X^^W) detg^'^'^'\X,a^,a2), 

i<*,i<^ 0=1,2 l<i<l 

l<j<r+s 

(88) 

with 

det ^(°'"''*)(A; ai; a^) = (aiaa)^^') JJ {X^~^ai - A'-^Os) I?("''^)(A), (89) 

l<i,jr<S 

and 

These facts are proved in all detail in Appendix |Al 

Before proceeding to the computation of the r.h.s. of eqs. fl87p . we come back for a 
moment to the argument we interrupted at the end of Section HI Eqs.f l87p have been 
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proven up to a constant independent of the difference n — k. To show that the constant 
is 1, it is enough to check the equations for t) and for i^2n+i(^j ^) hold true in the 

case n = k. For the l.h.s. we use eqs. (!56|) with q = e^'^*/^ and Zi = while for the r.h.s. 
we use eqs.dHMnO])- Rather than directly comparing the two sides of the equations, it is 
more convenient to compute the double ratio. A tedious but straightforward computation 
using Proposition 11041 shows that for both sides we have 



(91) 



which, combined with the direct verification for n = k = 1,2, gives the desired result. 
Proofs of the conjectures 

Taking the limit t — 1 directly in eqs. fl87l) is not easy. Instead we consider the ratios 
^g'^^fc "^hich are easier to compute, and from them recover eqs.( TT1l2ll4|6|) . Let us explain 

the computation for the case % i-t. ^ , the other case being dealt with in the same manner. 
Using the first of eqs. flSTl) we get 

g|„(fc-l,t) ,^u Y\T=t'-" t'"-^'-' - det (g(--"-^+^-^-^)({t^-n.^3^^^^2)) 

]-[2n-fc^2„,„fc_^i-l det(6;("'"-'='^)({t3-3|.^3^i^^2)) 

Then from eqs. fl88ti90|) we find that for generic values of v, A and Oj the ratio — detg( ■ • )(v,A,ai,a2) 



detg(^.'-+l'''-i)(v;A,ai,a2) 

does not depend on v and is given by a very simple formula 

det^{^,^+i>-i)(v;A,ai,a2) " ^ ' ""^^ V^^-^^^+^){\) ^^^^ 

J=-(s-l) 



p(-w)(A) ^ n:=i(A^-i)n-:7'(A^-i)n-;7'(A^-i) ^^^^ 

At this point we make use of these equations in eq. (l92|) and substitute \ = t^ , a\ = t and 
02 = t^. Repeating the same steps with the proper modifications for the other EFPs we 
finally obtain 

^2^(^-1;^) ^ ^..(n,.) (W\ [2n + fc-l],![n-fc],3![2A:-l],3![2fc-2],.! 
^n{^\t) V 3 ; [2n-fc],![n + fc-l],3![A;-l],3![3A;-2]i! 

g|n(fc-l;^) ^ ...-(n,.), ^ ^iSj^y"' [2n + - 1],! [n - fc],3 ! [2fc - 1],3 ! [2fc - 2],3 ! 

VSy [2n-fc]t![r2 + A;-l]i3![A;-l]i3![3A;-3]t![3A;-l]i 

£:2-+i(fc-l;t) _ ^3]* [2n + - k\,A\2k - \\^\\2k - 



S2n+i{k]t) V 3 y [2n-k + l]t\[n + k- l]t3\[k - l]ts\[3k - 2] J 

g2++i(fc-l;t) ^ /[3]*y~' [2n + fc],![n-fc + l],3![2fc-l],3![2fc-2],3! 

Sl+,{k-t) V 3 ; [2n-k + l]t\[n + k]t,\[k-iy\[3k-2],\ 



(95) 
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where we have introduced the usual t-numbers and t-factorials 



" — 1 

W- = and [i]t = j—j^ 

i=l 



The powers of t in the r.h.s. of eqs. fl95|) do not concern us because we are actually 
interested in the specialization t = 1, which at this point is immediate and reproduces 
the conjectured formulas ffH^lHlIB]) . 
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A A determinant evaluation 

In this appendix we evaluate the determinants of a family of matrices that appeared in 
the final step of the computation of the EFP in section m The matrices we are interested 
are labeled by three indices i, r, s and are made of blocks of rectangular matrices 



( <.(v) 












D 





(0) 



\ 



>2A) 



(96) 



where v = {wi, . 
lowing rectangular matrices 



vi), ttiX = {tti, ttiX, . . . , aiX^'~^^ ^} and each block consists of the fol- 



V 



V 



j+m-l j+m-1 



V 



j+m—1 



(97) 



The matrix ^(^'^''*)(v; A, ai, 02) has the total size 2{i + r + s) x 2{i + r + s). Here is an 
example 



?(2,1,1) 



(v; X,ai,a2] 





1 


1 


1 


1 











\ 




Vl 


V2 




Aai 
















vl 


vl 


al 


(Aai)2 
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1 


1 


1 


















V2 


a2 


Xa2 
















vl 


vl 


al 


(Aa2)2 
(Aa2)3 
(Aa2)^ J 




vl 


vi 


af 


(Aai)3 
(Aai)^ 


vf 


vl 


al 


I 


vf 


vl 


at 


vf 


vl 


al 
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We are interested in the determinant of Q^^'^'^\\^; \; ai, or equivalently in the deter- 
minant of the matrix ^(^'''''^^(v; A; ai; 02), obtained from ^(^'''''^^(v; A; ai, 02) through some 
simple row and columns manipulations 



^(^'^'^)(v;A,ai, 



02 J 



<1;.(V) <Wh-.(«iA) 



<l+2s;.H-.(«lA) i5ltH-2.;.(v) <U2s;.+. («2 A) 

(99) 

From the defining eq. fl99p we see that the first i columns of ^(^''"''^^(v; A, ai, 02) have rank 
i + r and therefore det ^(^•'■•''^(v; A, Oi, 02) = for r < 0. For r > 0, det ^(^''^■''^^(v; A, ai, 02) 
factorizes nicely. This property is the content of the following three propositions. 

Proposition 2. For r > we have 
dete?(^'^'^)(v;A,ai,a2) = J] (^'^ " ^i)' H (^* - A^" V) det ^?(°'^'^)(A, ai, 02) (100) 

l<ij<l 0=1,2 
l<i<i 
l<j<r+s 

Proof. The proof is done by induction on £. Let us look at det ^'^^''''^^(v; A; ai; 02) as a 
polynomial in vi. Using the matrix Q^^'^''^\-v; A; ai; 02) to compute the determinant, we see 
that its degree is at most 2(£+r+s— 1). The presence of all the factors [vi—Vj)"^, {vi — X^ai) 
and {vi — \^a2) is obvious using the matrix Q^^'^''^\-v; A, ai, 02). Since these factors exhaust 
the total degree it only remain to determine the term p(^''"'*)(v \ Wi; A; ai; 02) constant in 

Vl 

r+s-i e 
det^(''''^)(v;A;ai;a2) = J] {vi - X'ai){vi - X' a2)l[{vi - v,fV^'''■''\^r \vi; X; a,; a2). 

j=0 1=2 

(101) 

We can just set wi = in Q^^'^^'^^-v; X; ai] a2) and compute its determinant. We get 
immediately 

e 

detg^''^''\Y-X;ai;a2)U=o = X'^i^^lal+'a^+'Hv^ detg^'-^'^''\Y\vi; X;a,;a2). (102) 

1=2 

Comparing eq. fllOip for t>i = and ( I102p we get 

V'-'''''\y \ Vl- A; ai; 02) = det ^(^-^'^'^)(v \ v^- A; ai; 02) (103) 

which provides the recursion we were looking for. □ 

It remains only to establish the evaluation of det ^(°''''*)(A; ai; 02), but let us deal first 
with the particular case r = 0, which is quite simple and was useful in Section WA\ When 
r = the matrix ^^'^''''^^(A; ai; 02) has the form of a Vandermonde matrix and we find 
easily the following 

Proposition 3. 

det6;(°'°'^)(A;ai;a2) = (aia2)(^) J] (A*"' - A^"^)' J] (A'"'ai - A-^'^^aa) (104) 
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The case when r is generic is a bit more subtle and is dealt with in the following 
Proposition 4. 

det e?(°'^''*)(A; ai; aa) = (aiaa)^'^^ J] (A^-^ai - A^^^aa) V^'''\X) (105) 

i<«,i<s 



r'("''^)(A) 



(106) 

Proof. The proof is done by factor identification. For i = the matrix Q^^'^''^\X; ai, 
reads as follows 

Dl%{a^X) 



-Dr+2s;r+s(«lA) -D.^+2s;r+s(^2 A) 



(107) 



From the Laplace expansion of det ^'■'^■''''^-'(A; ai, 02) with respect to the first r + s columns 
we can easily deduce that all the determinants of the minors coming from the first r + s 

columns are divided by a} ^ , while all the determinants of the minors coming from the 

last r + s columns are divided by 03 • Iii this way we have obtained the factor {aia2p ^ ) 
in eq. fllOSp . In order to determine the total degree of det ai; 02) as a function of 

ai and 02 we notice that the only non vanishing contributions in the Laplace expansion 
of det ^''''''"'^•'(A; ai; 02) with respect to the first r + s columns come from the minor in the 
first r + s columns containing the first r rows and s among the last 2s rows, while the 
minor in the last r + s columns containing the rows from r + 1 to 2r and s among the last 
2s rows. Therefore the term of highest degree in ai corresponds to the minor containing 
the last s rows, and its degree is (''^'^) + s^. 

Now we show that as a function of ai, det ^^'^''"''^^(A; ai; 02) has zeros of order s — \j\ 
at ai = A-'a2 for < |j| < s — 1. Let for simplicity j > (the case j < being dealt 
with in the same manner). If we set ai = A-' 02 then the first r + s — j columns of 
-D^+2<i-r+s(A'''22A) are equal to the last r + s — j columns of D^^2s r+s{(^'^^) ■ Therefore if in 
^(°''^'")(A;ai = A^a2;a2) we subtract the r + s + j + i-th column from the z-th column, for 
l<z<r + s — j, we obtain a matrix whose first r + s — j columns have rank r. This 
means that the determinant of the original matrix has a zero of order at least s — j at 
ai = X^a2- 

Since we have determined all the zeros of det ^'^'^'^''^•'(A; ai; 02) as a function of ai and 
02 we have established eq. fllOSp up to the unknown factor T)^'^''^\\) which does not depend 
on ci and 02. In order to determine such a factor, we specialize ai = \^a2- We find that 
the first r columns of D^'^2s-r+s(A^a2A) are equal to the last r columns of Dfl2s-r+s{'^'2^) ■ 
Hence by subtracting the r + 2s + i-th column of the matrix ^'^"''''■'^(A; ai; 02) from its i-th. 
column, for 1 < z < r, we obtain the following matrix 



D^;,^(A^a2A) 

T, 

'r+2s\r+2s 



0P)(0) ( r \r+s \r+s+l \r+2s-l i \ X^+^-ll^l 



(108) 
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whose determinant is simply the product of the determinants of the two diagonal blocks 

det Dr°r and det -Dr+2s;r+2s 

det g?(0''-'^)(A; ai = X'a^; aa) = det D^^ det i5i%,.,+2s = 

By comparing eq. fllOQp with eq.f llOSp specialized at ai = A'*a2 we obtain eq. fll06p . □ 
Corollary 1. 

det g(^----)(v; A, ai,a2) . (.-i)(3.-2)/2 n ... . P^^-^^(A) ..... 

J— (-5 1) 



(111) 



^^^^-^•^^^HA) n:=i(A'-i)n-:7'(A^-i)n-;T'(A^-i) 

B Plane Partitions/Dimer coverings 

A Plane Partition can be seen as a tiling of a regular hexagon of side length with the 
following three types of rhombi of unit side length 









A /c-Punctured Cyclically Symmetric Self-Complementary Plane Partition (PCSSCPP) 
of size 2n is a Plane Partition symmetric under a 7r/3 rotation around the center of the 
hexagon of side length 2n and with a star shaped frozen region of size k (see Figure [1]). 

Following closely Ciucu [10] we compute the enumeration of /c-PCSSCPP of size 2n, 
that we call CSSCPP{2n,k). It is well known that Plane Partitions can be seen also 
as Dimer Coverings of an hexagonal graph. In the case of the /c-PCSSCPP the graph is 
reported in Figure [2j Thanks to the symmetry under a rotation of vr/S, it is sufficient to 
consider dimer coverings of the fundamental domain with "periodic boundary conditions" 
which are nothing else than dimer covering of a graph obtained by cutting the fundamental 
domain and joining the opposite cutted edges as on the right of Figure [21 Let us call 
G(n, k) this graph. Notice that G{n, k) is a planar graph symmetric under reflection along 
the vertical axis. In [10] Ciucu has proven that M{G), the enumeration of dimer coverings 



of a planar graph G with reflection symmetry, is related to the weighted enumerations of 
dimers covering of different graph G, obtained from G by removing the edges incident to 
the symmetry axis and lying on its right (see Figure [3]) • The relation reads 

M{G) = 2''M{G), 

where r is the number of edges lying on the symmetry axis, and the weighted enumeration 
M{G) is obtained by assigning a weight 1/2 to each dimer lying on the symmetry axis. 
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Figure 2: On the left the hexagonal graph corresponding to a PCSSCPP. The number 
of orange edges is 2?7, — /c , while the number of boundary vertices on the bottom side 
is 2n. The shaded region is a fundamental domain. On the right is the graph G(n, k) 
obtained by restricting to the fundamental domain and making the identifications of the 
orange edges imposed by the cyclic symmetry. 



In the case of PCSSCPP we are led to consider weighted dimers covering of the graph 
G{n, k) reported in Figure |31 which can also be seen as Rhombus Tilings of the domain 
D{n, k) or as Non-Intersecting Lattice paths starting from the right boundary of -D(n, k) 
and ending on its north-west boundary. This last representation allows to use the Lind- 
strom-Gessel-Viennot theorem and find 

CSSCPP{2n, k) = 2"-^M(G'(n, k)) = det[Q,+k,J+k]l<^,,<n-k (112) 

with 

The determinant of Qi+^j+k can be evaluated using Theorem 40 of [21] and we get 
CSSCPP(2n, k) = ff + 2^^1)1(3, +3^-2)1(3, +3^ ^2)! 



from which one easily finds 

CSSCPP{2n, k-l) {2k - 2)\{2k - l)\{2n + k - l)\{n - k)\ 



CSSCPP{2n, k) {k - l)\{3k - 3)\{3k - l){2n - k)\{n + k - 1] 



(115) 
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Figure 3: a) The domain G(n, k), obtained from G{n, k) by removing the edges incident 
to the symmetry axis and lying on its right. The dimers on the red edges, which lye on 
the symmetry axis, have weight 1/2. b) Another presentation of the graph G{n,k) and 
in dashed the corresponding domain D{n,k). c) The domain D{n,k). In dashed is the 
square lattice on which run the NILPs starting from the blue points and ending on the 
cyan points, d) A redrawing of the the square lattice of the NILPs. Each path going 
through a red edge gets a weight 1/2. 
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